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ABSTRACT 

We report the discovery of eight 7-ray pulsars in bhnd frequency searches using the 
Large Area Telescope (LAT), onboard the Fermi Gamma-ray Space Telescope. Five of 
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the eight pulsars are young (tc < 100 kyr), energetic {E > 10^^ erg s~^), and located 
within the Galactic plane < 3°). The remaining three are older, less energetic, and 
located off the plane. Five pulsars are associated with sources included in the Fermi- 
LAT bright 7-ray source list, but only one, PSR J1413— 6205, is clearly associated with 
an EGRET source. PSR J1023— 5746 has the smallest characteristic age (tc = 4.6 kyr) 
and is the most energetic {E = 1.1 x 10^'^ erg s""*^) of all 7-ray pulsars discovered so far 
in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages 
(tc ~ 1 Myr) and are the least energetic ~ 5 x 10^^ erg s~^) of the newly-discovered 
pulsars. We present the timing models, light curves, and detailed spectral parameters 
of the new pulsars. We used recent XMM observations to identify the counterpart of 
PSRJ2055+25 as XMMU J205549. 4+253959. In addition, publicly available archival 
Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023— 5746 
as a faint, highly absorbed source, CXOU J102302. 8-574606. The large X-ray absorp- 
tion indicates that this could be among the most distant 7-ray pulsars detected so far. 
PSR J1023— 5746 is positionally coincident with the TeV source HESS J1023— 575, lo- 
cated near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident 
with a 4.3(T excess reported by Milagro at a median energy of 35 TeV. Deep radio follow- 
up observations of the eight pulsars resulted in no detections of pulsations and upper 
limits comparable to the faintest known radio pulsars, indicating that these pulsars can 
be included among the growing population of radio-quiet pulsars in our Galaxy being 
uncovered by the LAT, and currently numbering more than 20. 



11 Subject headings: gamma rays: general - pulsars: general - pulsars: individual (PSR J1023- 

12 5746, PSR J2055+25) - X-rays: individual (CXOU J102302. 8-574606, XMMU J205549.4+253959) 

13 - open clusters and associations: individual (Westerlund 2) 



1. Introduction 



Until the launch of the Fermi Gamma-ray Space Telescope, only seven pulsars were found to 
have pulsed emission in 7 rays ( [Thompson 2008[ ), and only one of these pulsars (Geminga) was 
undetectable by radio telescopes (i.e. radio quiet). The Large Area Telescope (LAT) has detected 
7-ray emission from a large number of pulsars using radio ephemerides, both from young isolated 
pulsars (e.g. Abdo et al. 2009cpi ), and from millisecond pulsars (Abdo et al. 2009a). In addition, 
sixteen pulsars were found in blind frequency searches of the LAT 7-ray data ( Abdo et al.||2009"b ). 
Radio follow-up observations of those pulsars resulted in the detection of radio pulsations from 
^hree of them: PSRs J1741-2054 and J2032+4127 ( |Camilo et al.|2009D and PSR J1907+0602 ( [Abdo 



et al. 2010b). Two of these were found to have remarkably low radio luminosities. Perhaps more 



surprising is the fact that 13 out of these 16 pulsars were not detected in radio (Ray et al. 2010), 
and should therefore be considered radio quiet, or at least radio faint. 
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This paper reports on eight new 7-ray pulsars discovered in bhnd frequency searches of LAT 



ijiata, not included in the recently published First Fermi LAT Catalog of 7-ray Pulsars (Abdo 



et al. 2009d). As in the pulsar catalog paper, we present detailed timing and spectral results for 



each of these pulsars, including the pulse shape parameters, the fluxes, and the spectral indices 
and energy cutoffs for each pulsar. We also discuss possible associations, including previous 7-ray 



31 ijJetections. While a large number of the first sixteen 7-ray pulsars discovered in blind searches ( Abdo 



et al. 2009b) were found to be coincident with previously known 7-ray sources, such as EGRET 



33 unidentified sources, none of the pulsars in this sample have a definite counterpart in the 3rd 



34 EGRET catalog (Hartman et al. 1999), although one pulsar, PSR J1413— 6205, has a counterpart. 



EGR J1414-6244, in a revised catalog of EGRET sources ( [Casandjian k Grenier|2008[ ). This pulsar 



36 is also coincident with the AGILE source lAGL J1412— 6149 (Pittori et al. 2009). We also discuss 



39 



possible X-ray and TeV associations with the newly-discovered pulsars. The 8 pulsars identified 
here are among the 1451 sources included in the Fermi LAT First Source Catalog (Abdo et al. 
20T0a|). 



40 Finally, deep radio follow-up observations of these eight new pulsars have resulted in no new 

41 detections of pulsations, and we include the upper limits of our radio searches. 



2. Observations and Data Analysis 



43 The LAT is a high-energy 7-ray telescope onboard the Fermi satellite that is sensitive to 

44 photons with energies from 20 MeV to over 300 GeV. It features a solid-state silicon tracker, a 



cesium-iodide calorimeter and an anti-coincidence detector (Atwood et al. 2009). Events recorded 



with the LAT have time stamps derived from a GPS-synchronized clock on the Fermi satellite with 



< Ifis accuracy (Abdo et al. 20091). 



48 For pulsar science, the LAT represents a major advance over EGRET. It has a much larger 

49 effective area (~8000 cm^ for 1 GeV photons at normal incidence, or approximately 6 times that 

50 of EGRET), a larger field of view (~2.4 sr, or almost five times that of EGRET), and a finer point 

51 spread function (68% containment angle of 0?6 at 1 GeV for the Front section and about a factor 

52 of 2 larger for the Back section, vs. ~1?7 at 1 GeV for EGRET). The LAT also has a more efficient 

53 viewing strategy, operating primarily in continuous sky survey mode, as opposed to the inertial 

54 pointing mode used by EGRET This optimizes the amount of time the sky is in the field of view, 

55 covering the entire sky every two orbits (~3 hr). These improvements allow more photons to be 

56 accumulated per unit time, and provide a better signal-to-noise ratio in photon selection, making 

57 the LAT the first highly effective instrument for blind searches of 7-ray pulsars ( Abdo et al.|2009b ). 



^Available at 



http: //fermi .gsfc .nasa.gov/ssc/data/access/lat/lyr_catalog/ 



^For a single region in particular, however, a pointed observation would be more efficient for detecting pulsations, 
by accumulating more photons over a shorter time period. 
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58 For more details on the LAT, see Atwood et al. (2009) 



2.1. Blind Frequency Searches 

Since 7-ray photon data are extremely sparse (a typical 7-ray pulsar flux in the LAT may be 
of order 1000 photons per year), searches for 7-ray pulsars require long integration times. Further- 
more, young rotation-powered pulsars spin down as they radiate away energy, so their signals are 
not precisely periodic, making it necessary to search over a range of both frequencies and frequency 
derivatives. Fully coherent blind searches for 7-ray pulsars are therefore extremely computer inten- 
sive, since the number of frequency bins in the FFT increases with the length of the observational 
time (A^'bins = ST/max)) where T is the duration of the observation and /max is the maximum search 



frequency (Chandler et al. 2001). A fully coherent pulsar search, for example, would require the 



computation of hundreds of thousands of Gigapoint FFTs to cover the spin-down range of the ma- 
jority of young 7-ray pulsars. Such a search would also be highly sensitive to timing irregularities 
(e.g. timing noise and glitches). Instead, calculating the FFT of the arrival time differences, up 
to a maximum time difference of order ~week, greatly reduces the number of bins in the FFT. 
By doing this, we are able to reduce dramatically the computational demands (both in processor 
and memory), as well as the required number of trials to span the same parameter space, while at 
the same time reducing our sensitivity to timing noise. All of this is achieved with only a modest 



75 reduction in the sensitivity, relative to fully coherent techniques (Atwood et al. 2006). For these 



76 searches, we used a maximum time difference of 2^^ s (~ 6 days). We binned our difference search at 

77 a time resolution of 7.8 ms, resulting in a Nyquist frequency of 64 Hz for our searches and searched 

78 through 2000 steps of width A(///) = 5 x 10^^^ to cover roughly the spin-down range up to 

79 that of the Crab pulsar (-1.0 x lO'^^ s'^ < /// < Os"!). 

80 The data analyzed in this paper were collected from sky survey mode observations beginning on 

81 2008 August 4 (ME'I0= 239557414, MJD = 54682) and ending on 2009 July 4 (MET = 268411953, 

82 MJD = 55016). We used events with the most stringent background cuts (diffuse class photons of 



Pass 6.3 (Atwood et al. 2009)) with zenith angle < 105° (to avoid photons from the Earth's limb), 
and rocking angl^ < 40° . In addition, a few minutes were excised around two bright gamma-ray 
bursts (GRBs 080916C and 090510). We applied the time-differencing technique to photons with 
energies above 300 MeV, selected from a 0.8-0.9° (see Table 3) circular region of interest (ROI) 
around selected target positions given by the LAT first source catalog (Abdo et al. 2010a). We 



corrected the photon event times from each source to the solar system barycenter using the Fermi 



^Mission Elapsed Time (MET), the number of seconds since the reference time of January 1, 2001, at 0h:0m:0s 
in the Coordinated Universal Time (UTC) system, corresponding to a Modified Julian Date (MJD) of 51910 in the 
UTC system. 

*The rocking angle is the angle between the pointing direction of the LAT and the zenith direction, defined as the 
direction along a line from the center of the earth through the spacecraft. 
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Science Tool gtbary, assuming all photons come from the target position. We searched ~650 source 
positions from a preliminary version of the Fermi LAT First Source Catalog ( |Abdo et al.]|2010a| ) 
that were not coincident with likely active galactic nuclei (AGN). 

After finding a significant signal in the initial search, we followed up on the candidate signal 
by performing an epoch-folding search over a narrow region of frequency and frequency-derivative 



space using the PRESTO (Ransom 2001) pulsar software suite. All eight pulsars were discovered 
within ~ 11 months of routine survey observations. They have since been confirmed with several 
months of additional data. A pulsar detection is confirmed if adding the new data and folding using 
the original timing solution results in a continued increase in the significance of the detection of 
the pulsation, as measured by the chi-squared obtained from fitting the pulse profile to a constant. 

Compared to previous searches, the ones reported here are considerably more sensitive. In 
particular, these searches benefit from more than twice the amount of data, and a much better 
LAT source localization, than those that resulted in the discovery of the first 16 blind search 7-ray 
pulsars (Abdo et al. 2009b). The names and locations of the new pulsars are given in Table [!} 



103 jplong with known associations, including those from the Fermi LAT First Source Catalog ( Abdo| 



et al. 2010a) and the Bright Gamma-ray Source List (Abdo et al. 20091) 



2.2. Timing Analysis 

We have generated phase-connected pulse timing models for each of the pulsars. To do this, we 
first extracted photons for each source using a radius and energy cut chosen to optimize the signal 
to noise ratio for each pulsar. We corrected the photon event times to the geocenter using gtbary 
in its geocentric mode. We then determined a set of times of arrival (TO As) by first dividing our 
data into segments of approximately equal duration^ and then folding the photon times using a 
preliminary ephemeris to generate a set of pulse profiles. The TOAs were then measured by cross- 
correlating each pulse profile with a kernel density template that was derived from fitting the full 



113 mission data set, as described in Ray et al. (2010). Finally, we used Tempo2 (Hobbs et al. 2006) 



to fit the TOAs to a timing model that included position, frequency, and frequency derivative. In 
the case of PSR J1023-5746, our timing model also included second and third derivative terms, 
as well as a glitch, as described below. The timing solutions, along with the number of days of 
data and TOAs that went into generating them, and the rms of the resulting model, are given in 
Table [2] (with the positions in Table [T]). In the case of PSR J1413— 6205, there is an apparent glitch 
at ~ MJD 54735. Periodicity searches of the data before the glitch indicate that the magnitude 



^Note that at the time we carried out our searches, the 8 sources in which we found these pulsars were not known 
to contain a pulsar and their position was therefore much less well determined. 

®The number of days of LAT data used to generate each TOA varied for the different pulsars, ranging from ~9 
days for PSRJ2055+25 to ~23 days for PSR J1846+0919. See Table[2] 
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120 -was /^f / f = 1.7 X 10 ^, fairly large, though not unusually so for young Vela-like pulsars (Dodson 



et al. 2002 Saz Parkinson 2009). However, the short span of data before the glitch prevented us 



122 from including the glitch in the timing model and thus the model for J1413-6205 in the table is 

123 constructed only from post-glitch data. 

124 The case of PSR J1023-5746 is especially complicated. As is often the case in young pul- 



125 sars (Hobbs et al. 2004), this pulsar suffers from large timing irregularities. In addition to ex- 

126 periencing a glitch of magnitude A/// ~ 3.6 x 10^^ at ~ MJD 55041, this pulsar has a level 

127 of timing noise that required the use of higher order frequency derivative terms to model. Our 

128 model that fitted for position required terms up to / to obtain featureless residuals to the fit and 

129 resulted in a position that is ~2"from the Chandra X-ray source CXOU J102302. 8-574606. To test 

130 whether the proposed association with CXOU J102302. 8-574606 was compatible with the timing 

131 measurements, we made a fit with the position fixed at that of the Chandra source. This model 

132 required terms up to / , but also resulted in a good fit with essentially featureless residuals. The 

133 high order polynomial terms required to model the timing noise are strongly covariant with the 

134 position, resulting in the statistical errors being a significant underestimate of the true error. The 

135 Tempo2 timing models used in this paper will be made available online at the FSSC web sit^ 

136 Figure 1 shows the distribution in frequency and frequency derivative of the new pulsars (shown 

137 with unfilled triangles), compared with the previously known 7-ray selected pulsars (shown as solid 



138 triangles) (Abdo et al. 2009b) and the total known pulsar population. 



139 Figures [2jj9] show the folded light curves (two rotations, for clarity, and 32 bins per rotation) of 

140 the eight new pulsars in five energy bands (> 100 MeV, 100-300 MeV, 300-1000 MeV, > 1 GeV, and 

141 > 5 GeV), ordered by increasing right ascension. Each band illustrates the peak structure across 

142 different energies. The phases have been shifted so that the "first" peak lies approximately at 

143 0=0.25. The name of each peak (i.e. "first" vs "second") refers to their respective order of arrival 

144 relative to the off-peak emission. We caution, however, that for some pulsars (e.g. PSR J1954+2836) 

145 there may be no appreciable difference between the off-peak emission and the "bridge" emission 

146 between the two peaks, leading to some ambiguity in such designation. Consequently, there is also 

147 some ambiguity in our measurement of the corresponding peak separation. Five of the eight pulsars 

148 have a clear double peak structure, while three of them (PSR J1846+0919, PSR J1957+5033, and 

149 PSR J2055+25) are currently dominated by a single broad peak. For those pulsars showing a single 

150 peak, we tried folding the events at half or a third of the frequency (in case the original frequency 

151 found in our search was the second or third harmonic), but the various peaks obtained in this 

152 way were not statistically distinguishable. The additional statistics accumulated from continued 

153 observations by the LAT are necessary to firmly establish whether the single broad peaks are real, 

154 or whether they can be resolved into narrow multiple peaks with small peak separations. 

155 The pulse shape parameters, including peak multiplicity, 7-ray peak separation (computed by 

156 fitting the peaks with two gaussians and taking the difference between the means), and off-peak 



http : //fermi .gsfc .nasa.gov/ssc/data/access/lat/ephems/ 
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157 definition, are given in Table [3j 

158 2.3. Spectral Analysis 

159 We have performed the same spectral analysis on the eight pulsars as was performed on the 



160 46 pulsars presented in the First LAT Catalog of 7-ray Pulsars (Abdo et al. 2009d). The pulsar 



161 spectra were fitted with an exponentially cut off power law model of the form: 



KE^l^e ^cutoff (1) 



162 where the three parameters are the photon index F, the cutoff energy EcutoS, and a normal- 

163 ization factor K (in units of ph cm' -2 MeV"^), defined at an energy of 1 GeV. In order to extract 

164 the spectra down to 100 MeV, we must take into account all neighboring sources and the diffuse 

165 emission together with each pulsar. This was done using a 6-month source list that was generated in 



166 the same way as the Bright Source List, following the prescription described in Abdo et al. (2009d). 

167 We used all events in an ROI of 10° around each pulsar, and included all sources up to 17° into the 

168 model (sources outside the ROI can contribute at low energy). The spectral parameters for sources 

169 outside a 3° radius of the pulsar were frozen, taken from the all-sky analysis, while those for the 

170 pulsar and sources within 3° were left free for this analysis. In general, nearby 7-ray sources are 

171 modeled by a simple power law. In the case of PSR J1023-5746, however, the nearby sources include 

172 the 7-ray pulsar PSR J1028-5819 ( Abdo et al.||2009^ ), as well as the bright pulsar-like LAT source 



IFGL J1018. 6-5856, both of which were modeled by a power law with an exponential cutoff. For 
the eight pulsars, an exponentially cut off power law spectral model was significantly better than a 
simple power law. This can be seen by computing the test statistic TS'cutofT = 2Alog(likelihood) of 
the model with a cutoff relative to one with no cutoff, as shown in column 7 of Table [4j where the 
lowest value of TS'cutofT is ~20. Because the 7-ray sources we are fitting might contain a significant 
unpulsed component, we attempt to improve our fits to the pulsar spectrum by splitting the data 
into an on-peak and off-peak component, and performing the fit only to the on-peak events, while 
using the off-peak component to better estimate our background. Table [3] gives the definition of 
the off-peak phase intervals (the on-peak interval being the complement of the off-peak one), as 
determined by visual inspection of the light curves in all energy bands (the same method used in 
Abdo et al.||2009d ). We used the off-peak window to estimate the unpulsed emission. We re- fitted 



184 the on-peak emission to the exponentially cut off power law, with the off-peak emission (scaled 

185 to the on-peak phase interval) added to the model and fixed to the off-peak result. The resulting 

186 spectral index and energy cutoff for the on-peak emission of each pulsar are those given in columns 

187 4 and 5 of Table |4j For more details on the spectral fitting see Abdo et al. (2009d). 



Three of the pulsars (J1023-5746, J1413-6205, and J2055+25) showed significant emission 
{TS >25) in the off-peak component. The off-peak energy spectrum of PSR J1023-5746 shows no 
indication of a cutoff. Furthermore, the distribution of this emission does not follow the Galactic 
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plane (as might be expected if this were due to an improperly modeled diffuse background) and 
peaks close to the pulsar position, all of which suggests that such emission may be due to a pulsar 
wind nebula (PWN). However, the other two pulsars, J1413-6205 and J2055+25, have an off-peak 
component that is best fit by a power law with exponential cutoff, with photon index of F = 1 zt 0.6 
and r = 1±0.5, and £^cutofT = 1-8 ±1 GeV and -^cutoff = 0.65 ± 0.26 GeV respectively. This could be 



196 ipn indication of magnetospheric origin, similar to what is seen in the case of PSR J1836+5925 ( Abdo 



et al. 2009e). We caution, however, that these results are subject to some caveats. In the case of 



198 J1413-6205, the significance of the preference of the power law with a cutoff over a simple power 

199 law is marginal (a TS of 36 vs 26). The case of J2055+25, on the other hand, is highly sensitive 

200 to our definition of the off-peak phase interval. A reduction of the off-peak interval from 0.66 to 

201 0.5 (i.e. a 25% reduction) reduced the TS from 99 to 52 (nearly 50%), which seems to suggest 

202 that we may be overestimating the "true" size of the off-peak interval. After obtaining the spectral 

203 parameters for each pulsar, we compute the corresponding photon and energy fiuxes, and present 

204 the results in columns 2 and 3 of Table |4| It is interesting to note that the oldest pulsars tend to 

205 have the hardest spectra and lowest cutoff energies. Earlier studies in Abdo et al. (2009d) showed 

206 a weak trend of spectral index with E, and a stronger trend is evident here, albeit with lower 

207 statistics. 



208 3. The New 7-ray Pulsars 

209 The new pulsars are drawn from similar populations as the first 16 pulsars discovered in blind 

210 searches. Five of the pulsars (J1023-5746, J1044~5737, J1413-6205, J1429-5911, and J1954+2836) 

211 are energetic {E > 7 x 10^^ erg s~^) and young (tc < 100 kyr), with a large magnetic field at the 

212 light cylinder {Blc > 10 kG). The remaining three (J1846+0919, J1957+5033, J2055+25) have 

213 smaller Bic ~ 0.3-1 kG, are less energetic {E ~ (0.5 — 3) x 10^"^ erg s~^) and have much larger 

214 characteristic ages (tc > 300 kyr). 



3.1. Source Associations 



216 As shown in Table [T| five of the eight 7-ray pulsars are associated with sources found in 



217 the Fermi bright source list (Abdo et al. 2009f), while two of the pulsars were also reported as 



7-ray sources by the AGILE team (jPittori et al. 2009|). Only PSR J1413-6205 was clearly de- 

although it 



219 tected as a 7-ray source by EGRET (EGR J1414-6224 (Casandjian & Grenier 



2008 



is very likely that another one, PSR J1023-5746, was at least partly responsible for the EGRET 



221 source 3EG J1027-5817 (Abdo et al. 2009c). Three pulsars (PSRs J1023-5746, J1413-6205, and 



* Given the complicated nature of this region, it is safe to assume that at least part of the 7-ray emission from 
PSR J1413-6205 was captured as 3EG J1410-6147 ( [Hartman et al.|1999[ ), even if PSR J1413-6205 falls slightly outside 
the EGRET 95% statistical error error circle. 
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222 J1954+2836) may have even been detected as sources by the COS-B detector (2CG 284-00, 2CG 311- 



01, and 2CG 065+00 respectively (Swanenburg et al. 1981)), though it is difficult to establish a 
one-to-one correspondence between such sources and the pulsars, since the COS-B error radii are 
quite large (~ 1°) and source confusion becomes an issue. It appears, for example, that the COS-B 
source 2CG 284-00 was made up of various contributions which the LAT has now resolved into three 
separate 7-ray sources, including the pulsars PSR J1028-5819 and the newly-discovered PSR J1023- 
5746 dAbdo et al.|2009cD . We should point out that the issue of source confusion is one faced by the 
LAT too, albeit at a different sensitivity level than EGRET. The problem is particularly severe in 
the Galactic plane, where there is a large component of diffuse emission. In the case of PSR J1023- 
5746, for example, this resulted in an initial LAT source location that was far removed (> 10') from 



the true position of the pulsar (see Figure 11), making the early detection of pulsations from this 
source all the more challenging. 



234 3.2. Distance Estimates 

235 In the absence of distance information from observations in other wavebands, for example radio 

236 dispersion measures or optical parallax determinations, it is still possible to discern a general idea 

237 of the distance to these blind search pulsars. The method hinges upon the observed correlation 

238 between intrinsic 7-ray luminosity above 100 MeV and pulsar spin-down energy loss rate 

239 as depicted in Figure 6 of Abdo et al. (2009d). The luminosity trend is calibrated from the 



observed > 100 MeV fluxes using standard radio pulsar distance determination techniques and the 
presumption of a beam correction factoij^ = 1 for the 7-ray emission cone for all pulsars. 

The correlation is benchmarked at ~ 3.2 x 10'^^(-E'34)^/^ erg s^^ in Figure 6 of 



Abdo et al. 



(2009d), where £"34 is E in units of 10'^'^ erg s ^. 



Note that for E < 10^^ erg s"^ there is a break 



244 in the observed trend to a stronger dependence on E (e.g. L^ ~ (Abdo et al. 2009d). This is 



expected as pulsars cannot exceed 100% radiative efficiency. For the two pulsars that are below 
this transitional E (PSR J1957+5033 and PSR J2055+25), the efficiencies will be overestimated 
and thus the pseudo-distances given in Table |4] may be overestimates. 

There is considerable scatter (by factors of 3-10) of inferred pulsar 7-ray luminosities about 
this linear relationship, possibly due to the inaccuracy of the /n = 1 assumption. An earlier version 
of this correlation was calculated for the EGRET pulsar database ( Thompson et al.|1999 ). Physical 



251 origins for such a relationship are discussed in Zhang & Harding ( 2000 ) , having first been identified 



252 by Harding ( 1981 ) 



For a blind search pulsar with unknown dispersion measure, we can use this correlation, along 
with our determination of E from the measured spin-down parameters, to specify L^, and then 



^The beam correction factor, /n, allows us to convert the observed flux Fobs into 7-ray luminosity, L^,, using the 
formula L-y—ATrfnFobsD^, where D is the distance. 
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invert it to obtain a measure of its distance. This is referred to as the 7-ray pulsar pseudo-distance 
dps, and incorporates the presumption of /q ~ 1 in the absence of refined estimates of the sohd 
angle, which can be obtained using pulse profile information ( Watters et al.pOOQ ). In this way, one 
obtains 

dps = 1.6 kpc (2) 

Here Gioo-ii is the observed energy flux Gioo, as listed in Table [ij in units of 10~^^ erg cm~^ s~^. 
The significant scatter of inferred luminosities for radio-loud LAT pulsars about the correlation 
translates to uncertainties in pseudo-distances of the order of factors of 2-3. The resulting dps 
estimates for the new blind search pulsars are listed in Table 4, ranging from dps ~ 0.4 kpc for PSR 
J2055+2539 to dps ~ 2.4 kpc for the high E PSR J1023-5746. The relatively large pseudo-distance 
for PSR J1023-5746 is naturally expected due to interplay between its more-or-less typical Gioo 
flux and its high E. It is of considerable interest because of the possible association of this pulsar 



with the Westerlund 2 cluster, as discussed in Sections 4.2 and 14.31 



4. Multiwavelength Observations 
4.1. Radio 

These pulsars were all discovered in 7-ray searches and thus are 7-ray selected pulsars, but 
targeted radio observations are required to determine if they are also radio quiet, or could have been 
discovered in radio surveys independently. The population statistics of radio quiet vs. radio loud 



266 7-ray pulsars have important implications for 7-ray emission models (e.g. Gonthier et al. 2004). 

267 The precise positions derived from the LAT timing of these pulsars allowed us to perform deep 

268 follow up radio observations to search for pulsations from each of the new pulsars. We used the 

269 NRAO 100-m Green Bank Telescope (GBT), the Arecibo 305-m radio telescope, and the Parkes 

270 64-m radio telescope for these observations. The log of observations is shown in Table [5] and the 

271 instrument parameters used in the sensitivity calculations are shown in Table [6j All observations 

272 were taken in search mode and the data were reduced using standard pulsar analysis software, such 



as PRESTO ( Ransom|2001 ). We searched for pulsations for a range of values of dispersion measure 



(DM), from to some maximum, given in Table [sj 

To compute the minimum pulsar flux that could have been detected in these observations, we 



use the modified radiometer equation (e.g. Lorimer & Kramer 2004): 



,(S/iV)n,inTsys / W 



'-'mill — y „ I T — 7-, TiA v-'i 

Gy/upti^tAf V P-W 

275 where /3 is the instrument-dependent factor due to digitization and other effects; (S'/A^)inm = 5 is 

276 the threshold signal to noise for a pulsar to have been confidently detected; Tgyg = Tj-ec + ^sky; G 

277 is the telescope gain; Up is the number of polarizations used (2 in all cases); tint is the integration 



- 12 - 



278 time; A/ is the observation bandwidth; P is the pulsar period; and W is the pulse width (for 

279 uniformity, we assume W = O.IP). 

We use a simple approximation of a telescope beam response to adjust the flux sensitivity in 
cases where the pointing direction was offset from the true direction to the pulsar. This factor is 

/-(0/HWHM)2\ 

^=^^H 1.5 J' 

280 where 9 is the offset from the beam center and HWHM is the beam half-width at half maximum. 

281 A computed flux limit of S at the beam center is thus corrected to S/f for target offset from the 

282 pointing direction. The resultant flux limits are compiled in Table [5j We compare these flux limits 

283 with the measured fluxes of the population of pulsars in the ATNF pulsar catalog (fManchester et al. 



2005 ) in Figure 10 To make the fluxes comparable, we have scaled them all to the equivalent 1400 
MHz flux density using a typical pulsar spectral index of 1.6. Using the scaled flux densities and the 
pseudo-distances given in Table |4j we can estimate the luminosity limits at 1.4 GHz. Seven out of 
the eight pulsars have luminosity limits below L1.4 = 0.2mJykpc^. The remaining one, PSR J1023- 
5746, has a limit of ~ 0.6 mJy kpc^. Prior to the detection of radio pulsations from three LAT blind 



search pulsars ( Camilo et al.|2009 Abdo et al.|2010b ), the least luminous young radio pulsar known 



was PSR J0205+6449 in SNR 3C 58 with Li,4 « 0.5mJykpc2 ( [Camilo et aI1|2002[ ). Two of the 
new radio LAT pulsars, however, have luminosities an order of magnitude below this level. Thus, 
while the radio limits we present here are quite stringent by comparison with the overall known 
population of young pulsars, they are still far above the now least luminous pulsars known, and 
it is entirely possible that some of the new pulsars could be radio emitters at a level below our 
limits. It should also be noted that some of these pulsars, and in particular PSRs J1846+0919, 
J1957+5033, and J2055+25 located at high Galactic latitudes > 3°), may be very nearby and 



have very low values of DM. As was seen in the case of PSR J1741-2054 ( [Camilo et al. 2009), the 



received flux from such a nearby pulsar can vary greatly due to interstellar scintillation. The best 
way to detect such pulsars in radio might involve multiple low frequency observations. However, 
of the three pulsars mentioned, only PSRJ2055+25 has at present been observed more than once 
(see Table [5]). 

The lack of even a single radio detection among these eight pulsars brings the total number 
of known radio quiet pulsars to 22 (including Geminga), out of a total of 46 known young 7-ray 
pulsars, or ~50%. Such a high fraction may indicate that the size of the 7-ray beam is significantly 
larger than that of the radio beam, consistent with the predicitons of fan-beam outer-magnetosphere 
models for 7-rays and narrow polar-cap models for radio beams. However, more detailed population 
studies will be required to quantify this effect. 
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4.2. X-ray 



We searched archival X-ray observations for possible counterparts of the newly-discovered 
LAT pulsars, but found the pre-existing X-ray coverage of these fields to be extremely sparse. 
PSR J1023-5746 is the only pulsar with any significant X-ray observations, due to its proximity 
to the Westerlund 2 cluster. It was observed on three different occasion^^ (August 2003 and 
September 2006) by Chandra in the ACIS-S vf aint mode, with a total exposure of ~130 ks. In 
a new analysis of the Chandra data, we found a faint source at 02000 = 23™ 02;8, (52ooo = 
-57°46'07'.'01 with a 95% error radius of ~2", consistent with our best timing position for the 



316 pulsar obtained from LAT data (see Figure 12). This source is one of 468 X-ray sources previously 



317 reported by Tsujimoto et al. (2007) in their survey of the massive star- forming region RCW 49. 



318 While Tsujimoto et al. (2007) identify it as CXOU J102302. 84-574606. 9, it is also referred to in 



the SIMB AEF] astronomical database as CXOU J102302. 8-574606, and we hereon after adopt this 

320 designation. 

321 We find no optical source within 5" (upper limit my ~ 21) in the NOMAD optical cata- 

322 log. CXOU J102302. 8-574606 has a typical PWN-like power law spectrum with a photon index of 

323 1.2 ±0.2 and an Nh of (1.5±0.4) x 10^^ cm~^, resulting in an unabsorbed flux of 1.2ib0.3 x 10~^^ erg 

324 cm^^ s^^ (for a total of ~600 counts). No other simple spectral model was statistically suitable. 

325 The off-axis PSF of the Chandra ACIS camera is a complicated function of energy, and given that 

326 this source is ~ 8' from the center of the field of view we cannot make any claims about the spatial 

327 extension of the counterpart. The Nj^ value is comparable to the Galactic one (~ 1.3 x 10^^ cm~^) 

328 implying a fairly large distance (possibly greater than 10 kpc). The distance to the Westerlund 2 



329 cluster has been a matter of debate, with estimates as low as 2.8 kpc (Ascenso et al. 2007), or as 



high as 8.0 it 1.4 kpc (Rauw et al. 2007). The most recent estimate, based on an analysis of the 



CO emission and 21 cm absorption along the line of sight to the cluster places it at a distance of 6.0 
lb 1.0 kpc ( Dame|[2007 ). We note that the pseudo-distance that we have estimated for PSR J1023- 
5746 (given in Table 4) is 2.4 kpc. A rough estimate of the average transverse velocity of the pulsar 
from its birth site to its current location over its lifetime is given by v± ~ [(^8)('^2.4)/(''"4.6)]l-2 x 10^ 
km s^^, where 9s is the angular separation between the pulsar and its birth site, in units of 8 
arcminutes, ^2.4 is the distance to the pulsar, in units of 2.4 kpc, and t^ q is the age of the pulsar, 
in units of 4.6 kyr. This would be among the highest transverse velocities measured for any pul- 
sar (e.g. Hobbs et al.|2004 ). Accordingly, the association seems improbable unless a) the pulsar was 
born far from the cluster core (e.g. from a runaway progenitor), b) the cluster and pulsar distance 
are lower than our pseudo-distance estimate of 2.4 kpc, or c) the true pulsar age is substantially 
greater than its characteristic age. An on-axis Chandra observation might be able to resolve the 
expected bow shock structure of the apparent PWN and test this association. 



i°Obsids 3501, 6410, and 6411. 



http : / /simbad . u-strasbg . f r/ simbad/ 
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343 Shortly after identifying these 8 LAT sources as pulsars, we obtained short (~5 ksec) Swift- 

344 XRT observations in PC-mode of 5 of them (J1044-5737, J1413-6205, J1846+0919, J1957+5033, 

345 and J2055-I-25). In the first four cases, no likely X-ray counterparts were found within the LAT 

346 error circle. In the case of J2055-I-25, a 6.3 ks observation revealed two sources relatively close to 

347 the pulsar location (~l'away), but inconsistent with the current best position derived from the 



348 timing and likely associated with bright field stars (sources b) and d) in Figure 13). The upper 

349 limit on the flux for a putative X-ray counterpart in any of these five Swift XRT observations is 

350 1.5-2x10^^^ erg cm~^ s~^, under the hypothesis of a power law spectrum with a photon index 

351 of 2 and Nh of lO^i cm'^. Neither PSR J1429-5911 nor PSR J1954-F2836 have significant X-ray 

352 coverage. 

353 On October 26, 2009 the XMM-Newton satellite observed the field of PSRJ2055-F25 using 

354 the three European Photon Imaging Cameras (EPIC) MOSl, M0S2, and pn instruments. For 

355 MOSl and pn cameras thin filters were used while for M0S2 the medium filter was used. All 

356 three instruments obtained data in Full Frame mode which resulted in a time resolution of 2.6 s for 

357 the MOSl, M0S2 cameras and 73.4 ms for the pn camera. The data from each instrument were 

358 analyzed utilizing the XMM-Newton Science Analysis System software (SAS) version 9.0.0. The 

359 calibration data utilized by SAS was the latest available at the time of data reduction. We filtered 

360 the event data for bad events, retaining only those events from the MOS cameras with predefined 

361 patterns 0-12, and also excluded times of high background. 



362 In Figure 13 we show the combined image from the MOSl, M0S2, and pn cameras on 

363 XMM, resulting from an effective exposure of 19.0 ks after filtering. XMMU J205550. 8-^254048 

364 and XMMU J205547.2-F253906 (the "tadpole") are apparent in the Swift XRT image but owing 

365 to the higher effective area of XMM a number of fainter sources were also detected in the image 

366 including a faint source virtually coincident with the pulsar position: XMMU J205549. 4-1-253959. 

367 The region is crowded so in order to minimize contamination from other X-ray sources we form a 

368 spectrum by extracting events in a 29 arcsec circular region around the X-ray position centroid and 

369 also extract a background spectrum from a region that appears to be free of X-ray sources on the 

370 same CCD. For a circular region 29 arcsec in radius around a point source the MOSl PSF encircles 

371 roughly 85% of the photons at 1.5 keV. We can extract 97 photons from the source region and thus 

372 can form only an approximate spectrum for XMMU J205549. 4-1-253959. Fitting an absorbed power 

373 law to this spectrum we obtain an approximate fiux of ~ 2.4 x 10^^^ ergs cm~^ s~^ in the 0.5-10 

374 keV energy band with a power law index of F ~ 2.1. 



375 4.3. TeV 

376 Four of the eight new pulsars are in the Northern Hemisphere and four are in the South- 

377 ern Hemisphere. Only the Milagro TeV observatory, an all-sky survey instrument, has observed 

378 the four northern hemisphere pulsars. Milagro has reported a 4.3cr excess from the location of 



379 ^SRJ1954-F2836, with a measured flux at 35 TeV of 37.1±8.6 x 10"^^ TeV^^ s^^ cm"^ ( |Abdo 
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et al. 2009g). The remaining three pulsars in the Milagro field of view show no sign of TeV emis- 



sion, though this is perhaps not surprising given that they have spin-down luminosities that are 
30-200 times lower than PSR J1954-F2836. 

Of the Southern Hemisphere pulsars, three have no reported TeV observations. The remaining 
one, PSRJ1023-5746 is coincident with the TeV source HESS J1023-575 (see Figure [u]). The 
HESS collaboration reported the detection of TeV emission from a source with an extension of 
a = 0.18° lb 0.02 in the vicinity of the young stellar cluster Westerlund 2 (Aharonian et al. 2007). 



Three possible scenarios for the TeV emission were presented by the HESS team: the massive 
WR binary system WR20a, the young stellar cluster Westerlund 2, and cosmic rays accelerated 



389 in bubbles or at their termination shock and interacting with their environment (Aharonian et al 



2007). The estimated luminosity above 380 GeV (for an assumed distance of 8 kpc) is 1.5x10^^ 
erg s~^ or less than 1.5% of the spin-down luminosity of this pulsar, making the energetics of such 



an association plausible (e.g. Camilo et al. 2009). Recently, a jet and arc of molecular gas toward 
this source has been reported which suggests the possible occurrence of an anisotropic supernova 
explosion ( Fukui et al.|2009 ). The discovery of the very young and energetic pulsar PSR J1023-5746 
suggests that it likely plays an important role in the TeV emission from HESS J1023-575. Given 
the long list of known TeV PWNe (currently the most numerous class of TeV source^^ and the 
significant number of these associated with bright 7-ray pulsars, such a scenario seems at least 
plausible, if not probable. 



5. Conclusion 

The first five months of LAT data led to the discovery of the first 16 pulsars found from blind 



401 searches of 7-ray data (Abdo et al. 2009b). Here we report the discovery of an additional eight 

402 7-ray selected pulsars by performing blind frequency searches on 11 months of Fermi-LAT data. 

403 These new pulsars are largely drawn from the same population as the previous 16 (see Figure [s]). 

404 Deep radio follow-up observations of these newly-discovered pulsars suggest that they are all either 

405 radio-quiet, or extremely radio faint. 

406 While the less energetic 7-ray selected pulsars are found at high Galactic latitudes (|b| > 3°), 

407 pointing to a nearby Geminga-like population, the remaining ones, characterized by a higher E, are 

408 well aligned with the Galactic plane. Thus, the newly found pulsars represent well the two types 

409 of non-millisecond pulsars the LAT is unveiling: low E, nearby ones at high latitude, and high E, 

410 far away ones near the Galactic plane. Fainter, less energetic pulsars in the Galactic plane would 

411 be difficult to disentagle from the diffuse emission. 

412 Among the low latitude 7-ray pulsars, we find PSR J1023-5746, a neutron star with a remark- 

413 ably large rotational energy loss rate of 10^^ erg s~^, the largest so far of all 7-ray pulsars found 



^^http:/ /tevcat. uchicago.edu/ 
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414 in blind searches (Abdo et al. 2009b), and higher than >90% of the entire 7-ray pulsar popula- 

415 tion (Abdo et al. 2009d). The number of such energetic gamma-ray pulsars along the Galactic 

416 plane, however, could in fact be much higher. The very narrow distribution in Galactic latitude of 

417 the unassociated LAT sources (Abdo et al. 2010a) implies that many of them could be located at 

418 distances of several kpc, and have luminosities in excess of lO^^ergs"^. Using an average value of 

419 the conversion efficiency measured for the young 7-ray pulsars detected so far ( Abdo et al.||2009d ) , 

420 objects such as PSR J1023-5746 could account for many of the bright sources seen at low Galactic 

421 latitudes and so far unidentified. 
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Table 1. Names and locations of the new 7-ray pulsars 



PSR 


Source Association'^ 


R.A.'' 


Dec.'' 


r 


7 H 






hh:mm:ss.s 


dd:mm:ss.s 


( ) 


( ) 






1 0-93-09 9l'5'l 


-57-4fi-05('9'> 


984 9 


-0 4 




DFriT, 71094 0-^7^14 












HESS 71023-575 












r!X07T 7109309 8-574fi0fi 

K JW J \J ±\JjLikj\J^.O 1 '-t\J\J\J 












1Fr;7; 71044 5-5737 


1 0-44-39 R(^ ) 


-57-37-1 9 3('8'> 


986 6 


1.2 




1 Ar;7v 71043-5749 












1 Fr;7, 71413 4-fi205 


14-1 3-29 9l'1 "1 


-62-05-38C1 


312.4 


-0 7 




0FGLJ1413. 1-6203 












1AGLJ1412-6149 












EGRJ1414-6224 










J1429-5911 


1FGLJ1429. 9-5911 


14:29:58.6(1) 


-59:11:36.6(7) 


315.3 


1.3 




OFGL J1430.5-5918 










J1846+0919 


IFGL J1846.4+0919 


18:46:26.0(6) 


+09:19:46(11) 


40.7 


5.3 


J1954+2836 


IFGL J1954.3+2836 


19:54:19.15(4) 


+28:36:06(1) 


65.2 


0.4 




OFGL J1954.4+2838 










J1957+5033 


IFGL J1957.6+5033 


19:57:38.9(8) 


+50:33:18(9) 


84.6 


11.0 


J2055+25^ 


IFGL 2055.8+2539 


20:55:48.8(2) 


+25:40:02(3) 


70.7 


-12.5 



OFGL J2055.5+2540 
XMMU J205549.4+253959 



=^The current position uncertainty only allows for two decimal places in declination. 



''Sources are from the Fermi LAT First Source Catalog (IFGL, Abdo et al. 2010a), the 



Fermi LAT Bright Source List (OFGL, |Abdo et al.]|2009f| ), EGRET (EGR, [Casandjian 
Grenier|[2008l), and AGILE (lAGL, [Pittori et al]|2009p. We also list a TeV association 



(HESS) and X-ray counterparts identified with Chandra (CXOU) and XMM (XMMU). 

■^Right Ascension (R.A.) and Declination (Dec.) obtained from the timing model. The 
errors quoted are statistical (2.45 times the tempo2 uncertainties). They do not account for 
covariance between model parameters or systematic errors caused by timing noise, which can 
amount to several arcseconds, and should be considered when looking for counterparts. 

•^Galactic longitude (/) and latitude (5), rounded to the nearest decimal. 



-20- 



CO 



CO 

O 



< 

o 

Q 



CO 



i;D000001LOTt*CO!M 
COCN:— Hi-Hi-Hi-i:— ii-H 



lO O CN t~ CO 
Ol 00 00 t> 

o 



^ o o 
o 



CO 


CO 




C^j 








C32 




CO 




CO 


CO 


c3i 


I> 


CO 












CO 


CO 












CO 




00 




















CO 


00 


LO 


00 




CJl 


lO 




1— 1 


o 


o 


O 




CO 


o 


00 










I— 1 








CO 






I— I 


00 


CN 


o 


CO 






I— I 


CN 


I— I 


CN 




lO 


o 


LO 


I— ( 


LO 


CN| 


CO 


lO 


CO 


LO 




CO 


Ol 




00 




o 


LO 


CO 


CO 


CO 




CO 




LO 



CO CM lO ^ CO ^ 



^ 00 ^ CN ^ 
"J^ CO 00 CM lO CN O 
CT^ CTi CO lO 



OO 00 CN CN 
00 • • • • 
• CN (M (M O 

O 
CO 



o 
o 



(M O 



CN CN 



00 cm 

■CO UO 

1^ cm 

CN ^ 

00 

O CN 

cm 

cm ^ i-H 
00 cm 



o 

LO 

cm 

00 
CO 



00 



LO CM 

CO CO 

CO CO 

^ CM 



CO o 



LO 



CO 


1^ 


LO 


I— 1 




co 


o 


T— ( 




t~ 


CN 


im 


LO 


LO 


CO 


LO 


1 

CO 




1 

CO 


1 

<m 


CN 




1—1 


CM 


o 


o 






1—1 


1—1 


1—1 


1—1 





















co" 










o 










cm 


I— ( 


cm 


CO 


CM 


00 


CM 


LO 






00 




cm 


lO 




o 


CM 


CO 




00 


cm 


CO 


00 


CO 


CM 




1> 


CD 






CO 


CM 


CO 




1—1 






CM 


CO 




LO 




lO 




1—1 


o 


05 




t- 


1—1 


CM 






cm 


CO 


CO 






CO 


CO 




cm 


00 


o 


lO 


o 


CM 


LO 


CM 


+ 


+ 


+ 


+ 


CO 






LO 




LO 


LO 


lO 


00 


cm 


cm 


o 


1—1 


1—1 




CM 




1-5 


1-5 


1-5 



03 



CO 

LO 

00 

^ 

lO 

Q 

CO 

l>- 

LO 

I 

CO 
CM 

o 

1—1 

f-5 
Ph 



ft 

CD 
O 

X 



o 
o 

00 

LO 

Q 



CO 

s 
o 



CO 
1—1 

I— I 

1-5 

p^ 

CO 
PL, 



a; 



bjO 



CD 



O 



4J 
ft 

o; 

CD 

X 

CD 

co" 

1—1 

o 

LO 
LO 

CM 

00 
CO 

^ 

LO 

Q 

1-5 

o 

CD 
ft 

CD 



.5 ^ 



o 

bjO 



CJ 

o 
ft 

CD 

CD 
O 

a 

CD 

P-H 
=-t-H 



<1 



o 
o 
m 

LO 

Q 



bJO 



O 



LO 

Q 

1-5 

-d 

a 
j3 
o 



o 



o 

o 

CD 



o3 
O 



^ a> CO 



o 

LO 
LO 

I 

CO 

LO 



Q 

1-5 

CD 

bjO 

a 

Sh 

O 



LO 

CM 

+ 

LO 
LO 

o 

CM 
1-5 

P^ 

CO 

Plh 



LO 

o 

CO 







00 


CD 


_CJ 




CD 


(M 


Ph 




00 




1 


CD 


+H 


V 


requ 


O 
1— ( 

X 


03 


cteri 


p^ 
> 

CD 


-5746 : 


09(4) 


As th 


o3 

Ph 

03 

O 




CO 




O 


CD 




CM 


1 






O 


o 




+H 


o 










CO 


1-5 


II 


o 




A 


tPSR 




!-l 

0) 


•1— I 

bO 




nd 




00 




ive.' 


03 




X! 


las 


(M 
1 




an 


CD 




03 


CD 




CD 






t~- 


to 


> 


N 




ms 




CD 


K 


X! 




X) 




fl 


S 






1 


o3 




ed 


enc 


o 


del 


Co 


c/} 




X 


o 














CO 


CD 


cm" 




de 


CD 


X) 






o 


T3 


fl 




1 


a 




03 






he 


tanda 


>i 


II 


+H 


CJ 




CD 










<+H 


freque 


CD 


+H 


o 




03 


IB 

4H 




CD 


'3 




-(^ 




(-1 






CD 


bO 


CD 
fl 


'co 








CD 


CD 


■> 




1 


bC 


Ph 


ed \ 




+^ 
CD 


O 


bC 
fl 


n 






Xi 




obtai 


and 


del t 


ta use 


a 


CO 


CO 


o 




fl 


[» 


s 




!h 


O 






X! 


0) 


+^ 




o 






O 


a 




O 






. Colu 


atives 


gives 


ft 

o 


days 


CD 


X) 


> 




CO 




.2 




o 






CD 






nui 


CD 


X) 


fl 




ft 


>> 


spa 


olur 






O 


+^ 




s 


CD 


O 




O 


cr 


th 








CD 


CO 
CD 


o3 




rv 






CM 


CD 




> 


X! 




m 
Xi 


de 


'Sb 






O 


3 


LO 


CJ 




month 








e. Col 


third 


olumn 


:rom s 






X) 


U 




r nai 






X! 


;he 


o3 
X) 


CO 


CD 
o3 


03 

CO 




a 

o 


1 


Ph 

CD 


pul 


ove 


sec( 


Hz 


gen 



-S X! 

r^Sn CD 

bjO -o 



XJ 

03 
CD 

o 



^3 
O 

Ph 



-d 

03 



ft 

CD 



o 

Ph 
CD 
OJ 



CD 



bO 



^ ^ o 



CD 
fl 

bjO 
03 

a 

OJ 

> 

'bb 

o 



X! 

cm 

O] 

a 

o 
O 



Xi 

CD 



o 

Sh 



a 

o 



7^ ^. 



.-S 



o3 



bO 



X! O 

9^ y=l 



xi 



o 
xi 
I 

s 

•ft H 



CD 



Xi 

(1) 
bO 
03 



Ph 
CD 
Xi 

.a 

CD 



03 
CD 

CJ 



- 21 - 



Tabic 3. Pulse shape parameters for the new pulsars 



PSR 


Peak 


7-ray peak separation 


Reference phase 


off-peak definition 


ROI 




multiplicity 


A 


00 




(° ) 


J1023-5746 


2 


0.45 ±0.01 


0.66 


0.75-1.16 


0.8 


J1044-5737 


2 


0.35 ±0.01 


0.62 


0.66-1.16 


0.8 


J1413-6205 


2 


0.31 ± 0.02 


0.85 


0.69-1.09 


0.8 


J1429-5911 


2 


0.46 ±0.01 


0.86 


0.84-1.13 


0.8 


J1846+0919 


1 




0.36 


0.50-1.09 


0.8 


J1954+2836 


2 


0.43 ±0.01 


0.68 


0.84-1.16 


0.9 


J1957+5033 


1 




0.77 


0.53-1.13 


0.8 


J2055+25 


1 




0.83 


0.50-1.16 


0.8 



Note. — Light curve shape parameters of each pulsar, including the peak multiplicity, the phase 

separation between the two 7-ray peaks, the value of the phase at the reference epoch, barycentric 
MJD 54800 (UTC), and the off-peak region used in the spectral analyses. Column 6 lists the Region 
of Interest used to make the light curves, where we have selected only diffuse class events. 
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Table 6. Definition of observing codes^ 



Obs Code 


Telescope 


Gain 


Frcq 


A/ 




rip 


HWHM 


T 

rcc 






(K/Jy) 


(MHz) 


(MHz) 






(arcmin) 


(K) 


GBT-820 


GET 


2.0 


820 


200 


1.05 


2 


7.9 


29 


AO-327 


Arecibo 


11 


327 


50 


1.12 


2 


6.3 


116 


AO-Lwide 


Arecibo 


10 


1510 


300 


1.12 


2 


1.5 


27 


Parkcs-AFB 


Parkes 


0.735 


1374 


288 


1.25 


2 


7.0 


25 


Parkes-BPSR 


Parkes 


0.735 


1352 


340 


1.05 


2 


7.0 


25 


Parkes-lOcm 


Parkes 


0.67 


3078 


864 


1.25 


2 


3.3 


30 



^Tlie parameters in this table refer to Equation 3 in the text. 
''Instrument-dependent sensitivity degradation factor. 
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Frequency (Hz) 



Fig. 1. — The /-/ parameter space relevant to our blind searches. The pulsars in the ATNF cata- 
log ( Manchester et aL]|2005 ) are indicated by small black dots, and the radio-selected 7-ray pulsars 
are indicated by larger black circles (Abdo et al. 2009d). The 7-ray selected pulsars are indicated 



by triangles, where the filled-in ones denote previously reported blind search pulsars (Abdo et al. 



2009b), and the unfilled triangles denote the new pulsars reported in this paper. The shaded region 



corresponds to the part of the parameter space we have searched (see Section 2.1). The various 
dashed lines indicate how the different quantities derived from the timing parameters of the pulsars 
(characteristic age, magnetic field strength at the neutron star surface, and spin-down luminosity) 
vary across the parameter space. 



-26- 




Fig. 2. — Light curves for PSR J1023-5746. In this and the following figures, the folded light curves 
have a resolution of 32 phase bins per period. The black histogram in the second panel from the 
top shows the light curve for £^ > 5 GeV. The dashed line in the top panel represents the estimated 
background level, as derived from the model used in the spectral fitting. Two rotations are shown, 
for clarity, and the phase of the first peak has been placed at ~0.25 for clarity. 
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Fig. 3.— 



Light curves for PSR J1044-5737. 
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Fig. 4.— Light curves for PSR J1413-6205. 
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Fig. 5.— Light curves for PSR J1429-5911. 



-30- 




Fig. 6.— Light curves for PSR J1846+0919. 
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Fig. 7.— Light curves for PSR J1954+2836. 
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Fig. 8.— Light curves for PSR J1957+5033. 
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Fig. 9.— Light curves for PSR J2055+25. 
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Fig. 10. — Radio search upper limits (arrows), scaled to 1400 MHz using an assumed spectral index 



of 1.6, compared with the 1400 MHz flux densities from pulsars in the ATNF catalog (Manchester 



et al. 2005), shown as black dots. The three larger green circles represent the radio-detected LAT 



blind search pulsars, including PSRs J1741-2054 and J2032+4127 (Camilo et al. 2009), as wefl 



as J1907-I-0602, the pulsar with the lowest measured flux density, first detected in a 55 minute 



observation using the Arecibo 305 m radio telescope (Abdo et al. 2010b). 
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Fig. 11. — Fermi LAT image of a 1 square degree region of the sky around PSR J1023~5746. 
The smoothed counts map was generated using all diffuse class events above 100 MeV between 
4 August 2008 and 4 July 2009 (~ 21.2 Msec live time). The color scale is in counts per square 
arcminute. The overlaid green contours correspond to the significance {5a, 6.25cr, 7.5a, and 8.75a) 



reported by the HESS collaboration from the extended TeV source HESS J1023-575 (Aharonian 



et al. 2007). The LAT data are consistent with a point source. The white circle represents the 
i error circle) of OFGL J1024. 0-5754, the 7-ray source reported in the Fermi Bright 



position (and I 

Source List ( |Abdo et al.]|2009f| ). Note that the OFGL hst was generated using only three months of 
data. The blue cross (at R. A. =155. 806°, Dec. =-57. 743°, or ~2.1' away from the pulsar position) 
represents the preliminary estimate of the LAT source location based on 11 months of data and 
was the position used in the blind search that resulted in the discovery of pulsations. The large 
black rectangle is the area explored in the Chandra X-ray observations shown in Figure [T2j and the 
smaller black box within it, represents the 1 square arcminute region around the X-ray counterpart 



of PSRJ1023-5746, CXOU J102302.8-574606, also shown in Figure 12 at a larger scale. 
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Fig. 12. — Chandra ACIS-S (0.1-10 keV) X-ray image of the Westerlund 2 cluster region using 
~130 ks of data taken in August 2003 and September 2006. The color scale is in counts per square 
arcsecond. The small box on the right represents a 1 square arcminute region around the source 
CXOU J102302.8~574606, which we have identified as the X-ray counterpart of PSR J1023-5746. 
Note that the source is ~ 8' away from the core of the cluster. Inset - Zoomed- in image of the 
1 square arcminute region around CXOU J102302. 8-574606. The white ellipse represents the 95% 
confidence error ellipse of the position of PSR J1023-5746, based on the timing, as listed in TablejTj 
The errors are statistical only. Although the X-ray source appears to be extended, indicating the 
possible presence of a PWN, a full extended source analysis is dependent on a more complete 
understanding of the PSF of the instrument at such a large off-axis location. 



-37- 




Fig. 13.— XMM-Newton (0.2-12 keV) X-ray image of the region around PSR J2055+25, using ~26 
ks of data taken in October 2009. The color scale is in counts per square arcsecond. The green 
ellipse represents the current best location derived from pulsar timing, given in Table 1. North is 
to the top and east to the left. Four sources, labeled by letters, can be clearly identified in the im- 
age: a) XMMU J205548.0-h254115, b) XMMU J205550. 8-^254048, c) XMMU J205549.4-F253959, 
and d) XMMU J205547.2-h253906, the "tadpole". The best timing position of the pulsar is vir- 
tually coincident with that of XMMU J205549.4-F253959. XMMU J205550. 8-^254048 and XMMU 
J205547.2-I-253906 are apparent in a short Swift XRT image, but the source coincident with the 
pulsar, XMMU J205549.4-h253959, is not detected by the XRT. 



